Decoherence of adiabatically steered quantum systems 
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We study the efTect of Markovian environmental noise on the dynamics of a two-level quantum 
system which is steered adiabatically by an external driving field. We express the master equa- 
tion taking consistently into account all the contributions to the lowest non-vanishing order in the 
coupling to the Markovian environment. We study the master equation numerically and analyti- 
cally and we find that, in the adiabatic limit, a zero-temperature environment does not affect the 
ground state evolution. As a physical application, we discuss extensively how the environment 
affects Cooper pair pumping. The adiabatic ground state pumping appears to be robust against 
environmental noise. In fact, the relaxation due to the environment is required to avoid the accu- 
mulation of small errors from each pumping cycle. We show that neglecting the non-secular terms 
in the master equation leads to unphysical results, such as charge non-conservation. We discuss also 
a possible way to control the environmental noise in a realistic physical setup and its influence on 
the pumping process. 



I. 



INTRODUCTION 



During the past three decades along with a greatly 
increased interest in quantum mechanics and its appli- 
cations, we have witnessed a growing demand for accu- 
rate control of quantum systems. However, these systems 
are difficult to manipulate and very sensitive to external 
noise. 

A possible solution to the accurate control problem 
may be the adiabatic manipulation. This approach not 
only eliminates the need for fine tuning the time con- 
trol, but recently it has been proposed to solve com- 
plex computational problems^ and shown to be able to 
produce any unitary transformation.'^ Within the same 
framework, it is also possible to perform non-trivial peri- 
odic control to manipulate the quantum system by geo- 
metric means. This approach originates from the famous 
work of Berry'^ but later it has been extended to more 
general cases and analyzed for its possible applications 
in the so-called geometric quantum computing. ^^^ The 
main advantage is that the geometric dependence of the 
applied unitary operator renders it robust against fast 
fluctuations in the control signal. ^'^° The main drawback 
of this proposal is that to exploit all the advantages, the 
evolution time has to be long with respect to the time 
scale of the system. However, the evolution time cannot 
be increased indefinitely because the system undergoes 
dephasing and relaxation leading potentially to serious 
problems in maintaining the quantum state. Even if some 
efforts have been made to clarify the effect of the envi- 
ronment and its possible elimination in the context of 
geometric quantum computation, ^^'^^ a complete analy- 
sis is still missing. 

In this paper, we consider a two-level system initially 
in the ground state and study how it is influenced by 
the environment and an adiabatically changing exter- 
nal driving field. The time-dependent control makes the 



usual treatment^"^ by means of Markovian master equa- 
tion non-trivial.^^ Similar problems have been studied 
with different approaches"'^^""'^'' in the aim to estimate the 
decoherence effect induced by the environment. Here, we 
use a particular master equation derived taking consis- 
tently into account the drive, the environment and the 
combined effect. In particular, we include the usually ne- 
glected secular terms and find that, in our case, they be 
handled with extreme care since it can lead to unphysical 
results. Our result that zero-temperature environment 
does not affect the ground state adiabatic dynamics sug- 
gests that geometric quantum computation in the ground 
state^^ could be robust against decoherence. The analyt- 
ical derivation of the solutions of the master equation for 
different coupling and adiabatic evolution is discussed in 
detail. Similar results are obtained in the more realistic 
case of finite temperature environment. 

As a physical example, we study Cooper pair 
pumping^^ and discuss in detail the effect of the en- 
vironment in the pumping process. In this case, we 
find that the physical observables are robust during the 
ground state pumping. However, neglecting the non- 
secular terms in the master equation leads to a decrease 
in the predicted pumped charge and a difference in the 
average charge pumped through the two junctions, i.e., 
charge non-conservation. These results are obtained an- 
alytically by direct calculation and confirmed by numer- 
ical analysis. We propose a way to implement in situ 
the environmental noise by coupling a thermal resistor 
to the system through an array of SQUIDs. The main 
advantage of this proposal is that, by controlling the flux 
through the SQUID array, we can experimentally change 
the effective noise spectrum of the environment. 

This paper is organized as follows. In Sec. II the model 
is described. In Sec. HI, we write the master equation, 
solve it in the quasi-stationary limit, and discuss different 
approximations. In Sec. IV, we apply the present model 



to Cooper pair pumping in the presence of environment 
and discuss different regimes, the effect of temperature, 
and the imphcation of the secular approximation. Fur- 
thermore, we describe a way to engineer the environment 
in situ, and calculate the pumped charge. Section V con- 
cludes the paper. 



II. MODEL 

We consider a quantum system subject to external 
time-dependent control fields and interacting with the 
environment. The total Hamiltonian is H{t) = Hs{t) + 
He + V, where Hs{t) denotes the time-dependent sys- 
tem Hamiltonian with non-degenerate eigenvalues, He 
is the environment Hamiltonian and V is the system- 
environment coupling. We assume that V = Y (E) X, 
where X is the environment and Y is the system part 
of the coupling operator. The simplest way to study 
the dynamics of the system is to use the instantaneous 
eigenbasis of iJs(i), i.e., the adiabatic basis composed of 
states \tpn{t)) with energies En such that Hs{t) |^n(^)) ~ 
En{t) \ipn{t)) ■ We denote by D{t) the transformation 
from a given fixed basis to the adiabatic one. To prop- 
erly take into account the time-dependence of the basis, 
we study the dynamics of ptot = D\t)ptot{t)D{t) in the 
fixed basis. Here ptot{i) is the usual density operator of 
the total system in the Schrodinger picture. The evolu- 
tion of ptot is governed by the effective Hamiltonian 



ij(i)(t) = Hs{t) + hw{t) + V{t) + He 



(1) 



where Hs{t) = D^t)Hs{t)D{t), V{t) = D\t)VD{t) = 
Y{t) ® X, and w = -iD\t)b{t). 



We define the local adiabatic parameter as a{t) = 
h\\'w{t)\\/\/S.{t)\ where A(t) is the instantaneous min- 
imum energy gap in the spectrum and ||it;(t)|| = 
lTs\/w'^{t)w{t) is the trace norm oiw{t). This instanta- 
neous adiabatic parameter gives an accurate estimate of 
the local degree of adiabaticity of the evolution. Denoting 
with Tp the time interval for the adiabatic evolution and 
setting A„ii„ = min[A(i)], a rough estimate of a is given 
by a = h/ (AminTp) which represents a global adiabatic 
parameter. The w{t) scales as 1/Tp and then, for adia- 
batic evolution (a <C 1) it is usually neglected. ^^ How- 
ever, as discussed in Ref.^'^, it can be crucial to include 
its influence for the correct description of the dynamics 
of the system. In fact, the full Hamiltonian (1) suggests 
that the environment leads to relaxation to the instanta- 
neous eigenstates of Hs{t) + hw{t). If we neglect w(t), 
the environment produces relaxation in the adiabatic ba- 
sis of Hg{t) which usually differs from the correct one at 
the order 0(a). As we will discuss in detail, this contri- 
bution can be important in determining the expectation 
value of a physical observable. 

III. MASTER EQUATION 



Let us study the dynamics of a generic two-level system 
steered by a time-dependent Hamiltonian and subject to 
environmental noise. If the evolution is adiabatic and 
the system interacts weakly with the environment, we 
can treat hw{t) -f V{i) as a perturbation and derive a 
master equation for the density matrix in the interaction 
picture.^'* Up to the first order in a and the second in 
the system-environment coupling, we obtain 



dpijt) 
dt 



i[pi{t),wi{t)] ^2 







dt'[[pi{t)® pE,Vi{t%Vi{t)\ \ + 

[pi{t)®pEAwi{t')Mt")]]Mt)\ \, 



(2) 



r 



where pi{t) is the density matrix of the system in 
the interaction picture, Tr^ indicates trace over the 
environmental degrees of freedom and pe is the sta- 
tionary density operator of the environment. The in- 
teraction picture operators are defined as 0/(t) ~ 
e'"'^*/''Ulit,0)d{t)Us{t,0)e-'"''*/'', where C/s(t,0) = 
T'g^V' Jo Hs{T)dT jg ^]^g time-evolution operator, with T 
denoting time-ordering. In this frame, we can address 
and interpret the contributions of Eq. (2). The first con- 



tribution on the right is the driving term without system- 
environment interaction, the second one leads to the 
standard dissipative contribution of the Bloch-Redficld 
theory. The third contribution is a cross-term of the 
drive and dissipation. This last term is usually neglected 
as discussed above. ^^ 

Denoting the ground and excited state of Hs as \g) 
and |e), respectively, the master equation arising from 
Eq. (2) in the Schrodinger picture assumes the form 



Paa 



23m(w*gPge) + S{ujq) |m2| ^ - [S[-ujq) + ^(wo)] \m2\'^Pgg + 2 [5m (7712) 5m (pge) + 3?e (TO2) Sfte (pge)] 5(0)7711 



25(0) - 5(-wo) - 5(c^o) 



Wo 



{[3m (7772) Sm (TWge) + 3ffe (7712) 3fie (wge)] [Sm (7^2) Sm [pge) + 9?e (7712) 3ffe [pge)]} 



2S{Q) ^ S{-uoo) ~ S{u:o) 

Wo 

,5(0)-5(wo) 



7771 {3m (7712) 5m (wge) + 3fie (7^2) 3fie (wge)} Pgg 



Wo 



7n,i {3m (70,2) 5m (wge) + 5fte (7712) 5fte [wge)} , 



(3) 



and 



Pge = iWge{2pgg ^ l) + 7 ( Wge - Wgg)pge + iujQPge - S{L0a)mim2 + [^(-Wo) + 5'(wo)]m.i7772(Ogg - 2S{0)mlpge 

-^[^(-wo) + S'(a;o)]r772[3fie(7772)5m(/9ge) - 3m(7772)3fie(/9ge)] 



23(0) -S{-u;o)-S (loo) 



Wo 
5(0)-5(wo)^2 



Wo 



ra-^Wge — 777i2 



m^WgePgg 

S'(-wo) - 5(wo) 



Wo 



{3m(7(7,2)5fte(77;ge) — 3fie(7772)5m(77;ge)} 



25(0) - ^(-wo) - ^(wo) 

-2 777i{77n2[^m(Wge)Ke(pge) - Ke(Wge)^m(/9ge))J 

Wo 



I (7772) 3m(7i;ge) + 3fie (7772 )s3fie(wge)] Pge}, 



(4) 



r 



where /twq is the (instantaneous) energy gap between \g) 
and |e), 7771 = Ygg{t)/h = ~Yeeit)/h, 7772 = Yge(t)/h, and 
5(±wo) and 5(0) arc obtained from the power spectrum 
of the noise S{u}) = /_!^(X/(r)X/(0))e*"^dr.23 Here, the 
matrix element of a generic operator O is denoted by 
Oki = (fc|0|/) (with k,l = g,e) except for Wki ~ —i{k\l). 
Notice that the present equations generahze the ones pre- 
sented in Ref.^^ because they are obtained for a generic 
system-environment couphng operator. 

In the derivation of equations (3) and (4), we have 
imphcitly assumed that our system is in the Markovian 
regime. This imphcs that the bath correlation time r is 
short compared to the relaxation time of the system. In 
addition, since the evolution is adiabatic, wq, 7771, 7772, 
and w change slowly in time and we can approximate 

In writing the above equation we have neglected the 
corrections due to the Lamb shift. However, we expect it 
to have small influence on the dynamics of the system in 
the adiabatic limit. A detailed analysis is under way.^* 

Notice that Eqs. (3) and (4) arc, strictly speaking, 
valid only for adiabatic evolution and weak system- 
environment coupling. In addition, they are not in the 
standard Lindblad form and they do not a priori conserve 



the positivity of the density matrix. For this reason, par- 
ticular attention must be paid to the parameters used in 
the numerical simulations since they can lead to inaccu- 
rate results in the non-adiabatic case and strong coupling 
regime. 



A. Secular Approximation 

In Eqs. (3) and (4), we have included all the contribu- 
tions up to the second order in a and V contrary to the 
usually adopted secular approximation.^'^ In the interac- 
tion picture, the non-secular terms oscillate rapidly with 
respect to the relaxation-dephasing dynamics generated 
by the environment and thus their contribution averages 
to zero. Even though this approximation docs not lead 
to problems in many cases, we show that it is inadequate 
in describing adiabatic evolution. 

To compare the results with and without these terms, 
we explicitly write the master equations corresponding to 
Eqs. (3) and (4) but obtained with the secular approxi- 
mation as 



■sec 
rnn 



-25ni(w*g/9ge) + S{uja) jmsl ^ - [5'(-wo) + S{ujq)\ \m2\'^Pgg 

UJQ 

—2mi {5m (7712) 5m(wge) + iRe (TO2) SRe (w^e)} , 



(5) 



r 



and 



iWge{2pgg - l) + i{Wee - UJgg 

25(0) - S{-^o) - S{oJo) 



Wnn)pge + i^JoPge - [-S{-uJo) \m2\^ + 7:5'(a;o) |m2| ^ + 2S{0)ml]pge 



CJQ 



{2m2Weg + WgemDmiPge 

I 



(6) 



Notice that, in the time-independent basis (i.e. Wki = 0), 
these reduce to the usually adopted Bloch equations for 
a non-driven system 



aB 



Pgg - ^iPee - ^^Pgg, 

= iUjQPge — TgePge 



P% 



(7) 



where T 
and Ta 



I 



S{ll>o) |to2| 



t 



S{-ujo) \m2\ 



[^S{-iUo) \m2\ 2 + ^S{uJo) |m2| ^ + 2S{Q)ml 



These equations have immediate interpretation in terms 
of relaxation and dcphasing processes. 



B. Adiabatic quasi-stationary solutions for the 
ground state and the role of secular approximation 

In absence of system-environment interaction, i.e., 
S{uj) ~ 0, Eqs. (3) and (4) become the standard von 
Neumann equations for the density matrix 

Pgg == -23m(u;*gpge), 

Pge = iWge{2pgg -1) +i{Wee-Wgg)pge+ii^OPge-iS) 

Since we are interested in the evolution in the quasi- 
stationary limit for adiabatic evolution (a ^ 1), we look 
for the solution of pgg = and pge = 0. We have 



Pag = 
Pge = 



l + 0{a^), 



^ge 



+ 0(^2). 



(9) 



In general, Eqs. (3) and (4) can be integrated only nu- 
merically. However, we can obtain analytical information 
if we exploit the adiabatic limit. As above, the correction 
to the ground-state population pgg is of order a^ and the 
off-diagonal clement pge scales as a. 



In the following, we consider first the finite tempera- 
ture case, in which we assume that the excitation rates 
are exponentially small with respect to the relaxation 
ones; for example, Fge = ^eg cxp [— fia;o/(fcsT')], where T 
is the temperature of the bath. To the lowest order in 
a and in the quasi-stationary limit, the solution for pgg 
has a simple expression 



Pgg = 



'•g 



^ eg + t ge 



Oia^ 



0{a'). (10) 



The system remains in the ground state in the linear or- 
der in a, apart from leakage due to the finite temperature 
environment. 

Within the same approximation and using the above 
result, the equation for the off-diagonal term reads 

Pge = {iujo - 2S{0)ml)Age,i 



i^m {1112) 71128 {ujf))^e{Age^2) 
i^e [1112) m2S'(i:Jo)Sm(Age,2) = 0, 



(11) 



where we have defined 

(l - 23^^^^ ) w 
Age.l = — 



ffe 



CJo 



^Se,2 



1 + e" 



Pge, 

1 — e ^B^ 



Wo 



.(12) 



At low but finite temperatures, the solution of Eq. (11) 
is 






(13) 



Notice that, in the zero-temperature hmit, solutions in 
Eqs. (10) and (13) coincide with the the environment- 
free solutions (9): the ground-state evolution is robust 
and not influenced by the zero-temperature environment. 

This can be explained in a simple way.^^ The effective 
Hamiltonian in the absence of system-environment cou- 
pling is Hs + hw. In the adiabatic limit, the second term 
can be treated as a perturbation and the correction to the 
eigenstates of Hs can be calculated. This basis is usually 
called superadiabatic.^^ Up to the linear order in a, the 
superadiabatic ground state is \g') = {\g) — w* /wo|e)). 
In the adiabatic evolution, if the system starts in the 
ground state \g'{0)), we can assume that it remains in the 
eigenstate \g'{t)) and the corresponding density matrix is 
p' — l.9'(0)(5'(^)l- Thus, the density matrix elements in 
the Hs basis, pgg = {g\p'\g) and pge = {g\p'\e) satisfy 
Eq. (9). 

This result helps to understand the robustness of 
ground-state evolution subject to zero-temperature noise. 
In this regime, the environment cannot produce transi- 
tions towards the excited state |e'(t)) but induces only 
relaxation to the ground state \g'{t)). Since the system 
starts from the ground state and evolves adiabatically, 
it remains always in the ground state of Hs + hw and 
the environment has no influence on its dynamics. In 
this framework, it is clear that neglecting the w{t) terms 
leads to relaxation with respect to the basis of Hs and 
hence to different dynamics. 

The above result demonstrates the importance of the 
choice of the proper master equation. To emphasize this 
point further, we use a similar approach to solve the mas- 
ter equation arising from the secular approximation [Eqs. 
(5) and (6)]. In particular, in the zero-temperature limit 
and up to the order of the product of a and the square 
of the system-environment coupling, we obtain 



P] 



gg 

sec 



1, 



2iWn 



2iuJo — r ' 



(14) 



where T = S{ll!o) |m2| ^ -I- AS{0)'m1 and the corrections 
are of order of the product of the adiabatic parameters 
and the square of the system-environment coupling . The 
expression of the diagonal term is not influenced by the 
secular approximation, but the off-diagonal one depends 
on the coupling between the system and the environment. 
This is in contrast with both the physical picture dis- 
cussed above and analytical results in Eq. (13). This 
difference leads to different expectation values for physi- 
cal observablcs. 

We can obtain analytical information also in the op- 
posite regime in which the evolution is fast enough that 
non-adiabatic transitions determine the dynamics of the 
system. In this limit, neglecting the relaxation terms 
in Eqs. (3) and (4), we obtain again Eq. (8). The sys- 
tem is initially in the ground state with Pgg(O) — 1 and 
Pge(O) = 0. We expect that, due to sequential non- 
adiabatic excitations, in the quasi-stationary limit both 



the ground and the excited state arc equally populated, 
i.e., Pgg = pee = 1/2. With this assumption, the differ- 
ential equation for pge in (8) can be solved to obtain the 
quasi-stationary solution 



P99 = 1/2, 
Poe = 



(15) 



apart from corrections of the second order in the system- 
environment coupling. The system undergoes sequential 
transitions from the ground to the excited state and vice 
versa forming a fully mixed state. 

Thus, in the two limits in which relaxation dominates 
non-adiabatic transitions and non-adiabatic slight tran- 
sitions dominate relaxation, it is possible to find an ana- 
lytical quasi-stationary solution which correctly describes 
the system behavior [Eqs. (9) and (15), respectively]. No 
such solution is available in the intermediate regime in 
which relaxation and non-adiabatic transitions occur on 
the same timcscalo. 



IV. APPLICATION TO COOPER PAIR 
PUMPING 

In the previous section, we showed that ground-state 
evolution can be robust against environmental noise. The 
information about the dynamics of the system is con- 
tained in the off-diagonal term of the density matrix, see 
Eq. (13). Thus to test our theoretical model, it is nat- 
ural to focus on a physical observable which depends on 
Pge- An interesting candidate for the observable is the 
pumped charge in a superconducting circuit. Recently, 
it has been subject to intense study both theoretically^^ 
and experimentally^^ because of its connection to ge- 
ometric phases, ^^'^^'^^ and its potential application in 
quantum metrology. "'^^'^^^'^^ 

In particular, we consider the Cooper pair sluice^^ 
shown in Fig. 1. It consists of a single superconducting 
island, coupled to superconducting leads via two super- 
conducting quantum interference devices (SQUIDs). The 
SQUIDs operate as Josephson junctions whose critical 
currents can be tuned by magnetic fluxes. The electro- 
static potential on the island can be controlled by vary- 
ing a gate voltage, Vg, and there is a constant supercon- 
ducting phase difference, (p — ipn + ip^ between the two 
leads. The operator Uk — —idipi. (k ~ L, R) represents 
the Cooper pair number operator of the k-th SQUID. In 
the absence of noise, the Hamiltonian of the sluice can 
be expressed as^* 



Hs = Ec{n-ngf- Jl cos{ 



|-0)-JflCos(| + 0). (16) 



Here 9 = [pn. — '■Pl)/'^ and n = —ide are the operators 
for the superconducting phase and the number opera- 
tor of excess Cooper pairs on the island. The Josephson 
couplings to the left and right lead are denoted by Jl 
and Jn, respectively, Ug = CgVg/{2e) is the normalized 



(a) 




FIG. 1: (Color online) (a) Circuit scheme of the Cooper pair sluice (physical realization of an adiabatic two level system). The 
time-dependent control parameters are Jr, Jl, and Ug. The phase difference across the device is denoted by (p. The circuit 
in the dashed box is used to engineer the environmental noise characterized and modeled by a resistor R, a voltage 5Ve and 
coupled to the system with an effective capacitance Ce- (b) The time-dependence of the parameters {Jl, Ja, rit,} in a pumping 
cycle. 



gate charge, and Ec = 2e^/Cs is the charging energy of 
the sluice; Cg is the gate capacitance and Cs the total 
capacitance of the island. For further convenience, we de- 
fine the deviation of the gate charge from the degeneracy 



point as Sug 



1/2. 



The average-value given by the current through the 
k-th SQUID is Ik = — 2eTr(ptot'^fe) where the trace is 
over the degrees of freedom of both the system and the 
environment, and ptot is the total density matrix of the 
system and the environment. Using the von Neumann 
equation for ptot, we obtain 



2ie 

Ik = -T-JTr(ptot["fc,-ffs(i)]) +Tr(/9tot[nfe,ffi 

+Tript,t[nk,V]). 



(17) 



Writing the complete trace in terms of partial traces and 
tracing out only the environmental degrees of freedom, 
we obtain Tr5(p [uk, -ffs(t)]) which is the usual definition 
of the current for a closed system. ^^ The second term 
gives no contribution since Uk and He commute. The 
third term is the current directly induced by the environ- 
ment. Notice that, the last contribution is present only 
if [uk, V] y^ 0. In all the other cases, we can use the stan- 
dard definition of current operator: Ik = ^[nk,IIs{t)]. 
We study the case in which the noise is due to gate 
voltage fluctuations coupled capacitively to the sluice 
island. '^^ The charging part of the Hamiltonian in Eq. 
(16) has the form Ec[n - Cg/{2e)VgY. If an addi- 
tional noisy gate potential 5Ve is coupled to the sys- 
tem through a capacitance Ce, we can use the above 
expression to write the Hamiltonian describing the inter- 
action between the system and the environment. At the 
lowest order in 5Ve the Hamiltonian can be written as 
V = —2egn® 5VE{t), where g = CeJCy. is the coupling 
constant between the system and environment. Notice 
that since [''T.fe,'!^] = 0, the current operator of the fcth 
SQUID is determined only by the first term in Eq. (17). 



Its average value reads 

h = iPggIk,gg + Peeh^ee) + 2'?Re{pgeIk,eg) , (18) 

where Ik,nm = (nllkl'm) with m,n = g,e. The last 
term is the geometric contribution corresponding to the 
pumped charge 



Q^ 



2'^c{pgJk.eg)dt. 



(19) 



est charge 
can adopt 
denote the 
the island, 
sluice and 
where a^ 
identity. 



max{ Jl, J^j} and Ug w 1/2, only the two low- 
states are important for the dynamics and we 
the two-state approximation. Let |0) and |1) 
states with no and one excess Cooper pair on 
respectively. Thus the coupling between the 
charge noise has the form V = ega^ ® 5VE{i), 
= |0)(0| — |1)(1| a part terms proportional to 



A. Gate charge noise 

With the above notations and in this approximation, 
we can express all the relevant quantities in terms of the 
control parameters Jl, Jh and Ug as 



El2 

7 
V 

Wo 



J|, + 2 JlJr cos LP, 



1 
2 

arctan — tan — 

\Jr + Jl 2 



VK+d^)'' 



2E 



12 



h^T^' 



(20) 
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FIG. 2: (Color online) Pumped charge per cycle in a Cooper pair sluice for zero-temperature environment (a) as a function 
of phase ip with coupling constant g — Ce/Cy: = 0.01,0.0125,0.015,0.1 from bottom to top for / = l/Tp = 75 MHz, and 
(b) as function of coupling g with ip — 7r/2 and at frequencies / = 10, 25, 50, 65, 75 MHz from bottom to top in the region in 
which the arrow is plotted. In the numerical simulations, the parameters are J^' /Ec = 0.1, J^^ / jf' — 0.03 (with i = L,R), 
5nl' = -Sn"^ = 0.3, -Bc/fe = 1 K [Ec/{2Tih) = 21 GHz], R = 300 kfi, environment temperature T = 0, ^(ajo) = 2/iajo-R, 
S{~(jJo) = 0, and 5(0) = 2kBToR with To = 0.1 K. In inset we show the time evolution of the adiabatic parameter a{t) for 
/ = 75 MHz. 



and the instantaneous eigenstates assume the form 

Iff) = -L(yr^|o) + e-'^v/TT^|i)), 

|e) = -L(yri:^|0)-e-^v/T^|l»- (21) 

Further, the matrix elements of 1" {V — Y (E) X) read 
hrrii = ~grj and /im2 = g^/l — rf . 
The w matrix elements we need are 



Wn 



Wee 



Wn 



-2(l + '/)7 

I, 
-2(1 "'?)t 



(a/1 - '7^7 - i 



yr^^ 



(22) 



Using Ik = 2edip^Hs/h, the current operators are 
II = (2e/n)JLsin(f - 9) and Ir = {2e/h)jRSm{^ + 
6). We can write the current operators restricted to the 
{|0), |1)} basis using the formulas e*" = |1)(0| and e^** = 
|0)(1|. Inserting them in Eq. (19), the charge pumped 
through each junction reads 



n 



Qr 



2e 



dtJr 



ri'^c (pge) sin (^j + -j 



cos(^7+|j3m(pge)) 



dtJ, 



2e 
ft Jo 



'r]^c{pge)sm(^ 



■(cos (^7- -j3m(pge) 



^ 



(23) 



We define the average pumped charge as Q'~^ = {Q'f - 
Qp)/2 and charge conservation implies that Q'~^ = Q'£ - 

qI 



In our simulations the spectral density is taken to be 
ohmic such that S{uj) = 2hujR/[l - g-'^^/C^^^)] where R 
is the effective resistance of the environmental noise. For 
finite temperatures, we have S'(— w) = e~''"/''^'^"^-'5'(w) 
and S{0) = 2kBTQR, where Tq is the effective tempera- 
ture of the dephasing process. We note that the master 
equations presented in this paper are also valid for other 
types of environment spectra. 



B. Numerical Results 

The pumping cycle taken into consideration is shown 
in Fig. 1(b) and the corresponding evolution of the quan- 
tum system is obtained by numerical integration of Eqs. 
(3) and (4). The physical observable here is the pumped 
charge through the junction in the stationary regime 
which is reached after several sequential pumping cycles, 
i.e., after the initial transient is over. In fact, during the 
first cycles the pumped charge oscillates due to simul- 
taneous effects of non-adiabatic excitation and environ- 
mental relaxation and it takes several cycles to stabilize 
to the stationary values which are measured in the ex- 
periments. 

In Fig. 2(a), the pumped charge is shown as a function 
of the phase ip across the device ip and for different values 
of the coupling strength g while the pumping frequency 
is kept constant. The ideal pumped charge in the ab- 
sence of noisc"'^^ is also shown and it coincides with the 
numerical result in the strong coupling limit. For weak 
coupling to the environment, the pumped charge differs 
substantially from the ideal one but an increment of the 
system-environment coupling resumes the ideal pumping. 

In Fig. 2(b), we present the numerically simulated 
pumped charge as a function of g, for different pumping 
frequencies. For an adiabatic loop (e.g. / = 1/Tp — 10 
MHz), the pumped charge is close to the ideal value. This 
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FIG. 3: (Color online) Population of the ground state of J^s + fiui in the quasi-stationary evolution, (a) Evolution with g = 0.01 
and / = 10, 75, 100 MHz from top to bottom, (b) Evolution in presence of environmental noise for a non-adiabatic 75 MHz 
cycle with g — 0.01, 0.0125, 0.015 from bottom to top. The other parameters are the same as in Figure 2. 



confirms the first analytical results: in the adiabatic limit 
and T = 0, the environmental noise does not infiuence 
the evolution. For faster pumping frequencies, the be- 
havior is more complex because it is determined by the 
competing effects of non-adiabatic transitions and envi- 
ronmental relaxation. When the non-adiabatic transi- 
tions dominate over the relaxation process the predicted 
pumped charge is smaller than the ideal one. This is 
due to the fact that the excited state |e') of Hs + ^w 
is populated during the evolution and it pumps in the 
opposite direction with respect to the ground state |g'). 
Increasing the coupling we enter the regime in which the 
relaxation process dominates. Here, the system is forced 
to remain in the ground state \g') and we restore the 
ideal pumping. Thus, the numerical simulations confirm 
the analytical results obtained in Sec. Ill B for different 
value of a/g ratio. 

To check this interpretation of the numerical results, 
we calculate the population of \g') during the pump- 
ing cycles for different environmental coupling strengths. 
This is obtained by solving Eqs. (3) and (4) numeri- 
cally and projecting the obtained density matrix into 
the {IsO'k')} basis. Figure 3(a) shows how increasing 
the frequency leads to a leakage from the ground \g') to 
the excited |e') state. The figure shows that increasing 
the coupling strength from the region g ^ a <C 1 to 
a ^ g <C 1, the system begin to remain in the ground 
state l^'), thus restoring the evolution expected in the 
adiabatic limit without noise. This is a generic feature 
and we have chosen 75 MHz pumping frequency to illus- 
trate the phenomenon. 

Even if the convergence to the ideal evolution has 
strong physical motivation, we stress that the details of 
the numerical simulations in the extreme regime must 
be taken with care. In fact, in the intermediate regime 
both the analytical and numerical approaches fail to pre- 
dict the correct behavior of the pumped charge. In this 
regime, the numerical simulations in Fig. 2(b) are in- 
accurate and present an excess of pumped charge with 
respect to that in the adiabatic limit. This unphysical 
prediction is due to the fact that the master equation 



does not guarantee a priori the positivity of the density 
matrix and that we have truncated at the first order in a 
and g in the derivation of the master equation. Strictly 
speaking the g and a used are always small compared 
to 1 [see Fig. 2(b)] and thus, at least formally, they can 
be used as expansion parameters. Unfortunately, such 
analysis gives only an estimate of the corrections we are 
neglecting and proper analysis should include the higher 
order terms. ^^ 

However, there are indications that the features in the 
a and g regimes presented are realistic. The convergence 
to the ideal solution predicted by the analytical analysis 
is present in the simulations for all the range of param- 
eters used. It seems unlikely that such common behav- 
ior, which is realistic for small a and g, is a numerical 
artifact. The violation of positivity, which has a role 
in the intermediate regime, becomes less important for 
the maximum value of g presented. In fact, the analy- 
sis performed confirms that in these extreme regimes the 
positivity tends to be restored. 



C. Finite Temperature Pumping 



In Fig. 4, we compare Cooper pair pumping with zcro- 
and finite-temperature environments. We observe that 
the features of the zero-temperature simulation arc still 
present although, at finite-temperature, we obtain in gen- 
eral a smaller pumped charge with the same parame- 
ter values. This is expected since the presence of finite- 
temperature environment leads to an additional leakage 
from the ground state to the excited state, leading to a 
reduced pumped charge since the adiabatically pumped 
charge in the excited state is equal in magnitude but op- 
posite in sign to that in the ground state. ^^ 
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FIG. 4: (Color online) Pumped charge per cycle for finite-temperature environments (dashed curves) in comparison to the 
zero-temperature environment case (solid lines): (a) as a function of phase if for / = 75 MHz and g — 0.01, 0.0125, 0.015, 0.1 
from bottom to top and (b) as a function of the coupling strength g with (p = 7r/2 and frequency / — 10, 50, 75 MHz from 
bottom to top. The other parameters are as in Fig. 2 except for T = 0.03 K for the finite temperature simulations. 



D. Pumped charge with secular approximation and 
charge conservation 

All the results presented above depend on whether one 
includes in the master equation the non-secular terms 
and the combined effect of the environment and the driv- 
ing or not . From the analytical calculation presented in 
Sec. Ill A, we know that, in the quasi-stationary limit, 
the ofF-diagonal elements of density matrix within the 
secular approximation [Eq. (14)] have a different form 
compared with the full expression [Eq. (13)]. Thus we 
expect this to lead to different values of physical observ- 
ables. To compare the two methods, the pumped charge 
is again a privileged candidate since it depends directly 
on the off-diagonal matrix element pge- 

In the quasi-stationary limit, within the secular ap- 
proximation, at zero temperature, and for weak system- 
environment coupling, we obtain from Eq. (14) 



ge 



OJQ 



2a;o 



£2 

4a;, 



-—[^-^-^]=P^^-SP^ (24) 



where Sp denotes the correction with respect to the full 
solution. From this expression it results that p^g^^ is 
damped because of the interaction with the environment 
and this effect should reduce the pumped charge. To 
demonstrate this effect, we have chosen an asymmetric 
pumping loop in the flux and gate voltage parameters and 
performed a numerical integration of Eqs. (3), (4), (5), 
and (6). The results with and without the secular approx- 
imation are presented in Fig. 5. The pumped charge pre- 
dicted with the secular approximation decreases when the 
coupling increases whereas if we include the non-secular 



terms, we obtain a result which is essentially immune to 
g. Importantly, if the non-secular terms are included in 
the master equation, the two charges through the left 
and right junctions are equal but if we apply the secular 
approximation, we obtain a significant difference in these 
charges. The difference leads to charge non-conservation 
since no charge can be accumulated on the island for the 
quasi-stationary solution. 

We can also estimate analytically this charge non- 
conservation. Let us define the charge asymmetry be- 
tween the junctions as AQ ~ Ql'^'^'^ ~ Qr^'^'^j where 
the charge through junction k described 



-jG,scc 
-k 



by Eq. (14). Using Eqs. (20), (23), and Jr 



Jl sin (7 + I') / si^ (7 ~ ¥) ' ^'^ obtain 



AQ 



2e 



dt 



hfj 



- 2Ei2^m{Sp) 



(25) 



The first term in the integral gives no contribution since 
77(0) = ri(Tp). This is itself an important result since it 
means that, if we keep the non-secular terms (i.e. 6p = 
0), we obtain immediately that the charge is conserved. 
With the help of Eq. (24), the second term can be 
explicitly written as 



Qm((5/9) 



f^eiw. 



gel 



2^2 



Sm(wge)r2 



(26) 



Using Eqs. (20), (22), and the definition of F and S{uj) 
at zero-temperature given in Sec. IIIB and IV A, re- 
spectively, we can write explicitly the charge asymmetry 
as a function of time-dependent parameters 



10 



3/2 



AQ = 2eg-'h I dt[-i l-R {l - rf) + 



2Ei2 



^V9h{-R (1-,) ^ + 



^12 



(27) 
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FIG. 5: (Color online) Pumped charge obtained using differ- 
ent master equations presented for / = 10 MHz cycle. The 
straight line is the solution of the full master equation (3) and 
(4), the solid lines are the solutions of the master equation 
with secular approximation (5) and (6). In the latter case, 
the central solid line is the average pumped charge while the 
other two lines are the pumped charges through the left and 
right junctions. Here, the parameter cycle is strongly asym- 
metric with Jt/Ec = 0.1, Jl/Jl = 0.006, Jh^/Ec = 0.2, 
Jr/Jr = 0.04, (5nf = 0.4 5n™ = -0.03, T = 0, and the 
other parameters are as in Fig. 2. 



Ai = 



A2 



16{Ec5n'^Snyf 
m{Ec5n^5n'^'Y 

M{Ec5n^5nfY 

128(£;c(5n™5nf )5 

where AQ = Ei=i 5 Ai with Ay, = 0. The first two 
contributions are present only if the loop is asymmetric 



Aa 



A. 



(29) 



in the 5ng parameter (i.e. 



5nf 



^ liJ'T-^l) whereas ^4 



and ^5 yield non-vanishing contributions if the loop is 
asymmetric in the flux parameters (i.e. Jfj 7^ Jr)-^^ 
If the loop is asymmetric in both parameters as the one 
used in the simulations of Fig. 5, all the Ai terms in 
Eq. (29) contribute to AQ. 



E. Environment engineering 



In this expression, we can address five contributions. 
Since the third term depends only on 77 and its derivative, 
it yields no contribution when integrated along the cycle. 
The remaining terms can be calculated if we fix the phase 
to If — -K 12 and in the regime in which J™ ^ jf^ <C Ec 
and (Snf ,|<5n^'| > Jf^ jEc- Using Eqs. (20) and (22), 
we can write Eq. (27) in terms of the experimentally con- 
trolled parameters Ji, J^, and n^, and take advantage 
of the control cycle in Figure 1(b). To have an idea of 
how the final terms look like we write the first term in 
Eq. (27) as 



Ai = 



dt 



JrJl - JlJr) if {I - V^)^9^RSo 



2{Jl + Jly 



(28) 

It is convenient to divide the integration in six parts [see 
Fig. 1 (b)]: for example, in the first path only J^ depends 
on time. The remaining ?/ function can be approximated, 
since, for J™ ^ j/^ ^ Eq, it tends to a box function 
switching between ±1 [see Eq. (20)]. 

With this further approximation, we obtain four lead- 



The experimental set-up described in this paper has al- 
ready been used for the measurement of the Berry phase 
in adiabatic pumping of Cooper pairs. ^^'^^ To study the 
effects of the environment in the pumping, it is bene- 
ficial to control the coupling between the system and 
its environment. This can be achieved by adding to the 
already present intrinsic environment an experimentally 
controlled artificial noise source. Similar approaches have 
been developed in different physical systems for both 
testing the robustness of the studied quantum system 
and the cross-over between the quantum and classical 
behavior. ■^'^''^^ Here, we present a way to engineer in situ 
the environment noise in the superconducting circuits 
discussed above. 

The environment is controlled by the circuit in Fig.l 
which is composed of a thermal resistor R with noise 
voltage spectrum 5Ve- We place an array of m SQUIDs 
between the resistor and the island of the sluice as shown 
in Fig. 6. The array can be represented as a series of 
RLC parallel circuits with resistance Rs, inductance L, 
and capacitance Cs- The impedance of a single SQUID is 

y Rs Zcs ^L iL(tj>)u}Rs 

^RLC — RsZcs+ZlZcs+RsZl ~ iL{,p)uj+Rsil-L{<p)Lo'-'Cs) ' 

where L{(f>) = Lq/ cos{TT(t)/^o), Zl and Zcs are the 
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FIG. 6: (Color online) (a) Engineered environment connected to the sluice pump. The system-environment interaction is 
controlled by the flux through the SQUID array, (b) Pumped charge as a function of the flux through the SQUID array for 
Rs = 500 n, Cs = 0.3 fF (solid line) and Rs = I kn, Cs = 0.1 fF (dashed line). The other parameters are / = 10 MHz, 
LP = 7r/2, g = 0.025, m = 100 (number of SQUIDs in the array), Ci3 = 1 fF, Jc = 4 nA, i? = 1.5 kfl. In both cases, the 
simulations are performed with zero-temperature environment. 



impedances associated to the inductance L and capac- 
itance Cs, respectively, (/)/<I>o is the flux through the 
SQUID in units of the flux quantum, and Lq = h/(2TTeIc) 
where Ic is the maximum critical current of the SQUID. 
Note that the Joscphson inductance L{<f>) can be modu- 
lated by applying flux through the SQUID. 

The total series impedance of the circuit with respect 



J 



to the noise source in Fig. 6 is Ztot = R + Zqe +''nZj^Lc, 
and the voltage at the capacitor Ce is V = V ZcE/^tot- 
Here, we neglect the impedance of the sluice and pos- 
sible parasitic capacitance directly to ground inside the 
tunable environment. In this way, the gate voltage noise 
spectrum perceived by the system is 



~S{u:,cp) 



'-^Ce 



S{lu) 



S{u) 



L{(t>)uj'^R%CE-raL'^{<t>)uj*R%CsCE \ ' 



r 



_L2(0)a;2+_R^(l-L(0)a;2Cs)2 



RCeuj 



(30) 



Controlling the flux through the SQUID array, we can 
change the spectrum S{uj, </>) and hence we effectively 
tune the coupling between the system and the environ- 
ment. 

The pumped charge expected as a function of the flux 
through the SQUIDs in the array is presented in Fig. 6. 
The results are obtained solving the master Eqs. (3) 
and (4) numerically with the effective spectrum given 
by Eq. (30). This modulation of the the pumped charge 
could be large enough to be experimentally observable. 



V. CONCLUSIONS 

In conclusion, we have presented a detailed study of the 
effect of environmental noise on the adiabatic evolution 
of a two-level quantum system. In the zero-temperature 
limit, we find that the adiabatic evolution in the ground 
state is not influenced by the presence of the noise im- 
plying that the evolution is robust against environmental 



noise. 

These results are obtained taking consistently into ac- 
count the contributions in the master equation for the 
dynamics of the density matrix. In particular, we keep 
the non-secular terms and the terms describing the com- 
bined effect of the environment and the drive. This step 
is important to obtain the robustness of the ground state 
dynamics and, at the same time, to guarantee charge 
conservation. 

We have tested our theory by applying it to Cooper 
pair pumping in presence of system-environment cou- 
pling. The numerical simulations confirm the theo- 
retical prediction about the robustness of the ground- 
state pumping and suggest that in the non-adiabatic 
regime, relaxation can help to restore the desired ground- 
state evolution. In the same system, we studied the 
pumped charge in the secular approximation and ob- 
served that the pumped charge decreases with increasing 
system-environment coupling strength. Furthermore, the 
pumped charges through the first and the second June- 
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tion become different implying unphysical charge non- 
conservation. Thus we conclude that the secular approx- 
imation in this form cannot be pursued in adiabatic evo- 
lution. 

We have proposed a way to engineer the environmental 
noise in the discussed physical system. This allows us 
to modify the effective spectrum of the environment by 
changing the magnetic flux through an array of SQUIDs. 
The consequent change in the pumped charge could be 
experimentally observable. 

The results obtained concerning the robustness of 
adiabatic ground-state evolution are quite general and 
it would be interesting to extend our analysis. For 
example, adiabatic quantum computing^ is based on 
the adiabatic evolution in the ground state. Another 
possible extension is the construction of decoherence-free 
logical gates in the geometric quantum computation 
paradigm.^ In this case, the evolution occurs in a de- 
generate subspace and the manipulation of the quantum 
state is performed with the help of adiabatic evolution. 
In many of the proposals for geometric quantum gates, 
one of the main disadvantages is that the degenerate 
subspace is not the one with the lowest-energy. This 
likely leads to problems in controlling the relaxation 
to the effective ground state. However, geometric 
quantum computing can, in principle, be performed in 
the ground-state manifold. ^^ A detailed study of its 
robustness can be carried out by extending the theory 
presented here beyond the two-level approximation. 
However, this is left for future research. 



Appendix A: Derivation of the master equation 

The master equation (2) is obtained with a develop- 
ment taking into account all the terms up to the or- 



der wVV. As in standard derivation of the master 
equation^^, we assume that the density matrix of the 
environment is stationary and that the average of V over 
the environment degrees of freedom vanishes. 

Denoting the total density matrix of the system and 
environment as /5'°*(t) and employing the transforma- 
tion to the adiabatic basis as p\°^ = D^ p*'°^D, the von 
Neumann equation in the interaction picture reads 



pTHt) = j:[prit),hwi{t) + Vi{t)]. (Al) 

Notice that, in the weak coupling and in the adiabatic 
limit, wi{t) and V/(i) are perturbative contributions of 
different order. Tracing over the degree of freedom of the 
environment Eq. (Al) becomes 

/5Kt) = i[pi{t),wj{t)] + ^TrE{[pT\t),Viit)]}, (A2) 

where p^°* = Tr£;{p*°'}. 

Together with Eq. (Al), we employ the identity 



pr(r)=/5r(ri) 



dr'^rir'). 



(A3) 



Using iteratively Eqs. (Al) and (A3) we can obtain a 
perturbation expansion of Eq. (A2). 

Substituting Eq. (A3) in the last term in Eq. (A2) we 
have 

f>i{t) = t[pi{t),wi{t)] + ^TrE{[pT'{0),Vi{t)] 



Since the average of V over the environment degrees 
of freedom vanishes, Tr£;|[p^°'(0), V/(t)]} = and using 
Eq. (Al), we obtain 



J 



^i{t) = i[pj{t),wj{t)]-^TrE{ f dt'[[pr{t'),hwj{t%Vj{t)]} ~ ^Ti-e{ fdt'[[pr{t')Mt%Vi{t)]}-{A5) 



h^ 



h? 



The term pj°*(i') can be transformed using Eq. (A3); in the second and in the third term on the 
in particular, we substitute p\°*-{t') = P7°*(0) -I- right, respectively. Consistently, we have that 

J*'dt"^f\t") and pf\t') = pf\t) ~ jl,dt"f,f\t") Tr£{[[p*°'(0),w/(t')],V>(i)]} = and hence 



Mt) = ^[pJ[t),w,{t)]-^l:TE{ f dt'[[pr{t),Vi{t')],Vi{t)]} ^ ^TrE{ I dt' I dt"[[^r{t"),nwj{t')],Vj{t)]} 



ft nt 



h^ 



^^^e{J^ dt'J^^ dt"[[^r{nvj{t%vj{t)]}. 



(A6) 



Using Eq. (Al) we eliminate p\°^{t") from the above 
equation and, keeping the terms up to order wVV, we 



J 
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obtain 



t • t t' 

Pi{t) = i[p,(t),wiit)]-^TrE{l dt'[[pr{t),Vi{t%Vi{t)]} ~ ^Tte{ I dt' I dt"[[[pr{n,Vi{t'%wi{t')],Vi{t)\} 

+ ^TrB{ I* dt' 1^* dt"[[[pr{n, wj(t")],Vj{t')] , Vj{t)]}. (A7) 



r 



The third and fourth terms on the right are both of 
order wVV, namely, the highest order in our expansion. 
The last step in our derivation is to use Eq. (A3) to 

substitute p\°\t") with p\°\t) + j^" dt'" p\°\t"'); how- 
ever, since the terms with derivative of pj°* give contri- 
butions either of order w or V , they can be neglected. 
Thus, we can effectively substituted p\°^{t") with pj°'(t) 
in Eq. (A7) without introducing further approximations. 
The master equation (2) is obtained by rearranging the 



integration limits and the commutators of the last two 
terms of the resulting equation. 
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